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I. INTRODUCTION 
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We determine dynamical response functions {0\t,xi)0{'d,X2)) in the scaling limit of the quan- 
tum Ising chain on the half line in the presence of a boundary magnetic field. Using a spectral 
representation in terms of infinite volume form factors and a boundary state, we derive an expan- 
sion for the correlator that is found to be rapidly convergent as long as |2i±£2| > g 2 where ^ is 

■ the correlation length. At sufficiently late times we observe oscillatory behaviour of the correlations 
' arbitrarily far away from the boundary. We investigate the effects of the boundary bound state that 
' is present for a range of boundary magnetic fields. 
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. Over the last decade there has been significant progress in the calculation of dynamical correlation functions in 
I integrable quantum field theoriea^i^ii^i^i^i^. This has made it possible to determine the dynamic response of a variety 
of experimentally relevant one-dimensional models of Ising magnets^'-, integer— and spin-l/3ii quantum magnets, 
' I ' ■ Mott insulatorsi^, two-leg ladders^, carbon nanotube&i^ and ultra-cold atomic gasesi^. Having established the bulk 
D \ behaviour of such systems, it is of interest to investigate the effects of impurities on the dynamic response. The 
^ simplest possible effect of non-magnetic impurities in a spin chain material is to break the chains into finite segments. 

Resulting changes in the susceptibility have been investigated in both spin-l/2l^ and spin-1^^ Heisenberg chains as well 
I t . as two-leg Heisenberg laddersi^. In the latter two cases the chain break can lead to the formation of edge states, which 

■ have been observed in inelastic neutron scattering experimentsi^. Impurity and boundary effects are also relevant 
^ , to scanning tunneling microscopy experiments^, where the associated breaking of translational symmetry can lead 

' to fingerprints in the dynamic response that assist in characterising of the underlying bulk behaviour—. For the 
aforementioned reasons it is important to extend the calculation of dynamical response functions in one-dimensional 
models to systems with boundaries. In the gapless case this is readily accomplished by boundary conformal field 
theor y^^i^^ , while the gapful case is as usual more involved. One successful approach has been the extension of the 
truncated conformal space approachSi to systems with boundar y^^i^^ . The method we will follow in the present 
work is the extension of the form factor bootstrap approach to boundary integrable field theories2ii2^i22ii2£i2i. We will 
focus on the simplest integrable field theory, the quantum Ising chain- i^^i'^^i^'^i'^^i^^i'^^i^^ . One-point functions in the 
• boundary Ising model have been studied by Konik et al.?^ . The purpose of the present work is to analyse two-point 
I functions relevant to scattering experiments. 
00 • This article is organised as follows. First, we will review the basic facts on the Ising field theory in the bulk. We 
' proceed by discussing the Ising model on the half line in the presence of a boundary magnetic field and the form factor 
: bootstrap approach to boundary integrable field theories. In order to set the stage, we calculate the Green's functions 
' of the Majorana fermions, which show typical light-cone effects. After a brief discussion of the local magnetisation we 
C , proceed with our main result, the calculation of the dynamical spin-spin correlations and the corresponding spectral 
' function. The derived expansion in terms of infinite volume form factors and a boundary state is found to be rapidly 
J> , convergent as long as |^i±^| > 0.2 where ^ is the correlation length. As for the Green's functions we observe light-cone 
effects, which result in oscillatory behaviour of the correlations arbitrarily far away from the boundary. Furthermore, 
^ . , we study the effects of the boundary bound state, which exists for sufficiently small values of the boundary magnetic 
' field, on the correlation function. Finally, we discuss the two-point function of the disorder field, which shows the 
same qualitative features as the spin-spin correlations. 



II. QUANTUM ISING CHAIN 



The Hamiltonian of one-dimensional quantum Ising model is given by 

i^latt = -J^(fT?a?+l+A<). (1) 

i 

Here a'^ and cr^ are the Pauli matrices and J is the exchange energy. The model ([T]) is related to the two-dimensional 
classical Ising model, see for example Ref. [3§. The critical temperature of the latter is related to the coupling constant 
of the quantum Ising chain by A = 2^^^. The Hamiltonian ^ is invariant under the Z2-transformation af — > af, 



(T? — * —erf - For A < 1, which corresponds to the low-temperature phase of the 2D classical Ising model, this symmetry 
is broken. The order parameter field takes a non-zero expectation value and the ground state is two-fold degenerate. 
For A > 1 the model has a unique ground state. This regime corresponds to the disordered high-temperature phase 
of the classical 2D Ising model. The point A = 1 is the location of a quantum phase transition. 

At small deviations from criticality, |A — 1| ^ 1, one can pass to the continuum limi t^^''*" . This leads to the Ising 
field theory, which is defined by the Euclidean action^i^ 

S = — I (fx {ijdtp + i)dij -t- i M{pi}) , (2) 
27r J 

where ip and 7p are the two components of a Majorana fermion. The complex coordinates used in ([2]) are given by 
z = T + ix and z ~ r — ix, where t = it denotes imaginary time. This results in d — dz = (dr — idx)/'2 and 
d = dz ^ {dr + idx)/'2. The mass is proportional to the distance from the critical point, M ex. J(l — A). The model 
(12) is conformally invariant at the critical point M = (see for example Ref. l4Th . In the ordered phase, which we 
will consider throughout this paper, the mass is positive. Furthermore, we set the velocity to one and use the short 
distance normalisations 

z(0|V'(T,a;)V'(0,0) |0) ^ 1, Izl^/"^ (0| ct=^(t, a;) cr^(0, 0) |0) ^ 1, as \z\ ^ 0. (3) 

For later use we will need the notion of mutual semi-locality of operatorsii^i^i^^. Let us consider the operator 
product Oi(t, x) 02(0, 0) = Oi{z, z) 02(0, 0). If we take Oi counterclockwise around O2 in the plane, i.e., we perform 
the analytic continuation z — s- e^'^' z, z — > e^^'^' z, the operators Oi and O2 are said to be mutually semi-local if 

Oi(e2"'z,e-2-z)O2(0,0) =/oiO.Oi(z,z)O2(0,0). (4) 

The phase I01O2 is called the semi-locality factor. The two fields are mutually local if I01O2 — 1- Semi-locality is the 
mildest form of non-locality, in general the right-hand side of ([4]) may be more complicated. The mutual semi-locality 
factor of the spin and disorder operators can be read off from their operator product expansion^i^ 

a^{z, z) f,^{0, 0) ^ {e'-/'V~zm + e-'^/^V^zm)- (5) 

This implies that when taking once around /i^ one obtains an extra minus sign, i.e., l^^^^ = —1. In the same way 
one finds Z^^z = l^^^ — l^a-^ — l^^^ — ~1. On the other hand, the disorder field is local with respect to itself. 

In order proceed further, we need to construct a basis of scattering states. In order to do so, we need to specify 
which "fundamental field" (that is a field that has a non-vanishing matrix element between the vacuum and the 
one-particle states) we will use to create the fundamental excitations. In order to avoid additional signs in many of 
the equations it is customary to choose the fundamental field to be bosonic. We will therefore use the disorder field 
fi^. This implies that the fundamental excitations are viewed as bosons with two-particle scattering matrix = — 1. 
Let us denote the corresponding annihilation and creation operators by A{9) and A^{9) respectively. They fulfil the 
Faddeev-Zamolodchikov algebra 

A{e,)A{02) = sA{e2)A{0,), 
A^iei)AHe2) - sA\e2)A^0i), (6) 

A{ei)AH92) = 2TTS{0i~02) + SA\e2)A{9i), 
where the scattering matrix is 5 = — 1. The vacuum state is then defined by 

A{9)\0) = 0, (7) 

and a basis of scattering states is given as 

\ei,...,e,-,)=AH9,)...A^9,,)\0). (8) 
In terms of the states (H)) the resolution of the identity reads 

id^\o){o\ + Y^^ / ^^-^\e,,...,9^){e^,...,e,\. (9) 
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For the calculation of correlation functions the knowledge of the matrix elements or form factors of local operators 
is necessary. In the ordered phase the non-vanishing form factors of contain an even number of particles and are 
explicitly given by^'^^'^'^ 

2" - 

/(0i,...,02„) = (Ok^(O,O)|0i,...,02„) =iVo J^tanh^y^, (10) 

where cto = (Olf^ |0). For comparison, the form factors of the disorder operator [i^ are only non- vanishing if the 
number of particles is odd, 

2n+l fl _ 

(O|a^^(O,O)|0i,...,02„+i) =i'Vo n tanh^^. (11) 

i<j 

The form factors (jlOp and pTjl follow from the following set of requirementsi'^'^: 

1. The form factors /(^i, . . . ,6*71) are meromorphic functions of 9^ in the physical strip < 3m 0„ < 27r. There 
exist only simple poles in this strip. 

2. Scattering axiom: 

/(^i, . . . , 6i+i,6i, . . . , On) — S f{9i, . . . , 9i,9i+i, . . . , 6,1), 
with the scattering matrix S* = — 1. 

3. Periodicity axiom: 

fiOl +27:1, 02,..., 0n) = lo^^ f{02, ...,en, 9l), 

where O = or . The mutual non-locality phases are given by Zcr>^ = — 1 and /^^^^ = 1, respectively. 

4. Lorentz invariance: 

/((?i+a,...,0„+a) = e^"/(ei,...,0„), 
where s denotes the spin of the fields, which is s^r^ = and s^^ = as well as = —1/2 and = 1/2. 

5. Annihilation pole axiom: 

Res[f{e',B,ei,...,en),e' = e + iTT\ =if{e^,...,e,,) [i-^o;.^nL^ • 

Note that the squared brackets on right-hand side equal 2 for as well as /i^, as the extra minus sign due to 
la^fi'^ = — 1 is compensated by an additional factor S. As there exist no bound states in the Ising model, these 
are the only poles of the form factors. 

We note that if we had used the Majorana fermion tp as fundamental field, the excitations would be viewed as fermions 
with unity scattering matrix. Furthermore, there would appear additional minus signs in the axioms above^^. 



III. BOUNDARY ISING MODEL 



We now turn to the Ising field theory on the half-plane (r, x), t e M, x G (— oo,0]. The boundary is located at 
X — Q and r denotes imaginary time (r = it). The Hilbert space of states associated with the semi- infinite line 
r = const., —00 < X < 0, is denoted by Tie- The boundary condition is imposed through application of an external 
time- independent magnetic field h which couples to the boundary spin crg(T) = t7^(T, 0). The action for this system 
is given by22 

5=2^/"dT/" dx{i:dij + '4:di; + iMij^l;) + ^ 1 dr fhal-^-Hj-^adro) . (12) 
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A^Oi) A^{e2) A^e) 

FIG. 1: Two-particle scattering and scattering off the boundary. 



It was shown by Affleck and Ludwig^'? that in the critical model the application a magnetic field at the boundary 
generates a flow from free to fixed boundary conditions. The properties of the boundary spin operator ct^ were 
analysed in Refs [23ll43l . In the Lagrangian framework defined by p2)) it can be written as21 
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'^bir) ^ ^iHr,x)+^PiT,x)) air) 



x=0 



(13) 



Here a(r) is an additional fermionic boundary degree of freedom which is introduced to describe the ground-state 
degeneracy and anticommutes with ip and ip. The cases h = and h oo correspond to free and fixed boundary 
conditions, respectively. We note that the application of a non-zero boundary magnetic field removes the ground-state 
degeneracy of the Ising model. 

In terms of the Majorana fermions the boundary condition reads 



i-(^-v;) 



/j2 

x—0 A 



(14) 



As was shown by Ghoshal and Zamolodchiko\«22 the Ising model with boundary conditions (|14p still possesses infinitely 
many integrals of motion (which can be constructed from the integrals of motion of the model in the bulk) and hence 
remains integrable. 

If one thinks of the boundary as an infinitely heavy impenetrable particle B sitting at a; = and writes |0b) = 
S |0) G Tie for the ground state in the presence of the boundary, the scattering of particles created by the operators 
Al' {&) at the boundary is encoded in the relations 



A\e)B = R{e)A\-B)B. 



(15) 



The function is interpreted as the single-particle reflection amplitude off the boundary. In order to preserve the 
integrability of the corresponding bulk system, the boundary scattering matrix R{&) has to satisfy several conditions 
as discussed in Ref. [13, which leads to the following result for the boundary scattering matrix for the Ising model 
with boundary magnetic fleld h 



R{0) 



, , ITT U 

1 tanh 

4 2 



K — i sinh 
K + i sinh 9 ' 



K = 1 - 



2M' 



(16) 



Free boundary conditions are recovered for h — whereas fixed boundary conditions are obtained in the limit h oo 
(we restrict ourselves to ft, > 0). The purpose of the present work is to calculate the two-point correlation function 



CiT,xi,X2) = (0B|r,a^(r,a;i)a^(0,X2)|0B), 
where T-r is the time-ordering operator. The time-dependence of the operators is given by 

a'{T,x)=e^"^a'{0,x) e""^^. 



(17) 



(18) 



where Hb is the Hamiltonian of the system in the presence of the boundary. 

As we are working in the Euclidean formalism, r and x are interchangeable. Therefore, one may also take x to 
be the Euclidean time, which implies that the equal-time section is the infinite line, x ~ const., — oo < r < cxd, 
and the associated Hilbert space H is that of the corresponding bulk theory The boundary at a; = now 

appears as an initial condition which is expressed in terms of a "boundary state" |B). As was shown by Ghoshal and 
Zamolodchikov^^, the correlation function (fT7|) can then be expressed as 



C(t, Xi,X2) 



\%a^{T,xi)a%0,X2)\B) 
(01 B) 



(19) 
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a%T,Xi) 
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FIG. 2: Euclidean rotation. The boundary condition at s = turns into an initial condition. The "time" in the rotated system 
runs between a; = and a; = — oo. 




FIG. 3; K{^) as function of ^. Note that h = 2vM corresponds to «: = —1 and h = he = v2M to «: = 0. 



Here 7^ is the ai-ordering operator, which orders the largest Xi to the right, and 
model on the infinite line. The boundary state is given by 



where K{£^) = R(i-k/2 — £). For the Ising model it is explicitly given by 



G 7i is the ground state of the 



(20) 



This simplifies to 



i^te\ ■ ^ 1 K + COsh^ 



iffroc(C) = -i coth |, infixed (C) = i tanh |, 



(21) 



(22) 



for free and fixed boundary conditions, respectively. For later use we introduce the real function K[£^) — — iif(^), 
which is plotted for several values of k in Fig. [3l We note that for < /i < 2y/M the function K{S,) is negative for all 
^, whereas for 2\fM < h the function K{(,) is positive in a certain ^-interval. Furthermore, for < h < he — V2M 
one finds K(^) < —1 for certain ^. In this region the boundary scattering matrix ()16|) has a pole in the strip 
< 3m < tt/2, indicating the existence of a boundary bound state^^. We will discuss this state in more detail in 
Sec, nil A1 below. Finally, we note that for free boundary conditions the boundary scattering matrix (|16p has a pole 
at ^ = 0. This results in the appearance of a zero-momentum mode in the boundary state, i.e., pop is replaced by^l 



|Bf,ee) = (1 + A\0)) exp ( ^KiOA^^OAHO ] |0) 



(23) 



Since all form factors of involving an odd number of particles vanish, this zero-momentum mode does not contribute 
to the two-point function of cr^, and hence will be ignored in our analysis. Finally, we note that K(£^) — —K{—£^) 
and K{£^) — > (^ — > 0) except for free boundary conditions. 



5 



FIG. 4: Contour of integration for the excited boundary state. 



As we have interchanged space and time and x is running from to — cx3 in the new framework, the r- and 
x-dependence of operators (r, x) is now given by 



a"(r, x) = e~^" e"'"^ (7^(0, 0) e'"^ e"^, 
where H is the Hamiltonian of the system on the infinite hne, — oo < t < cx3, and P is the total momentum. 



(24) 



A. Boundary bound state 

If we consider free boundary conditions {h = 0) there exist two degenerate ground states |0b,±). The effect of 
a small field < h < he = y/2M is to split these into two non-degenerate states |0b) and |1b), where |1b) can be 
interpreted as a boundary bound stated. In this domain wc can parametrise k as 



K = cosv, < V < 



(25) 



The two states can be distinguished by the asymptotic behaviour of the one-point function (OB|cr^(a;) |0b) +<yo 
or (1b| cr^(a;) |1b) —fro for x — oo, respectively (we assume that h > 0). The energy of |1b) is given by 
El — Eq + Msinw, where Eq denotes the ground-state energy. If h approaches the critical value he, the boundary 
bound state becomes weakly bound and its effective size diverges. For h > he ^6]) possesses no pole in the physical 
strip and hence no boundary bound state occurs. 

Correlation functions in the boundary bound state can be calculated following Ref. [13 



(IbI r,a^(r,a;i) (7=^(0, X2)|1b) = 



{0'\%a^{T,xi)a-^{0,X2)\B') 



(0'|B'> 

where |0') = |0, — ) G 7i denotes the "wrong" ground state. The excited boundary state is given by^I 



IB') 



exp 



(26) 



(27) 



where the contour of integration is shown in Fig. [H The contour encircles the pole of (j2ip at ^ i u, whose residue 
equals — 2i cot u tan(w/2). 

Due to the Z2-invariance of the Ising model in the bulk the matrix elements of remain unchanged up to a minus 
sign when considering the state |0'), i.e., the form factors are 



(0 {^1 , ■ • ■ , ^2n }') = {0'\a^A\ei)...A^ (02„ ) 1 0') - - 



2n ^ 

i"cro Y\_ tanh- 

i<j 



(28) 



In particular, we have Eq^ = E^,. The minus sign can be deduced from (0'| a" |0') = - (0| a" |0). 
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a) 



b) 



oscillating behaviour 



FIG. 5: Light-cone eflFects observed for the Green's function pip , a) In the bulk one observes oscillating behaviour if < . 
b) The boundary contributions show oscillating behaviour if 4R^ < . 



IV. GREEN'S FUNCTIONS 



In order to set the stage, we first calculate the Green's functions of the Majorana fermions. The mode expansions 
in the presence of the boundary and in the absence of a boundary bound state are given by 



OO 

M r de 



2 J V27r 





OO 

M f de 



i Mx sinh 6 



-i Mx sinh 8 



M X sinh 8 



2 J V27r 





c{e)e 



T cosh 9 ^ 



gf giA/a:sinhe ^ ^j-^-j ^-f ^-iMisinhe 



gf g-iMa:sinh6» 



I g-f gi A/2:sinhe 



Here c{9) and c^(6') are canonical fermion annihilation and creation operators {c(0), c^{9'y\ — 2ttS{9 ~ 9'). 
The Green's function of ■0 is now easily derived to be (we assume r > for simplicity) 



(OB|V'(T,xi)0(O,a;2)|OE 



M 

T 

M 



d9t 



g-e gi A/rsinhe 



— / d9R{9 + i9o)e- 

^ .1 — , — , 



(29) 



(30) 



(31) 



MV4:R^+t'^ coshe 



where = arctan(2|_R|/r) and Ki denotes the modified Bessel function of order one''^. Furthermore, we have 
introduced centre-of-mass coordinates R — {xi + X2)/2 and r — X2 — xi. We stress that R < 0. In real space, t = it, 
the first term is oscillating for < and damped for > t^, i.e. we observe a light-cone effect. The second term 
is oscillating for AR^ < and damped otherwise. The physical interpretation of the oscillating behaviour is that for 
4i?2 < & particle can propagate from (0, a;2) to {t,xi) via the boundary (see Fig. [S]). If one calculates the Green's 
function using the rotated system and the boundary state, one obtains an additional phase of 7r/2. The physical origin 
of this phase is the Lorentz spin = —1/2 of ij:, which implies that the Green's function transforms non-trivial under 
Lorentz rotations. In the case h < there is an additional term due to the presence of the boundary bound state. 
In the same way one obtains (we assume r > 0) 



(Ob I V'(t, xi) ^5(0, X2) |0b) = -i A/ ifo (Mv/;;^T^) - i ^ 

z \ IT — Zri 



d9R{9 + i9Q) e 



-MVUF+t' cosh 6 



(32) 



which shows the same kind of oscillatory behaviour as ([HI 
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FIG. 6: Local magnetisation in the ground state and the boundary bound state for different values of the boundary magnetic 
field > up to first order in the boundary reflection matrix K. The full lines represent the ground-state result p3p . the 
dashed lines the result for the boundary bound state (|34|> . 



V. LOCAL MAGNETISATION 

In this section we briefly discuss the one-point function of the spin operator. The one-point function in the ground 
state I Ob) was first calculated by Konik, LeClair and Mussardc^? and to second order in K is given by 



(OBk^(a;)|OB) -ao-iao^"^X(0 tanh ^ e^*^^ ^^i^ « 
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^^I'^l nif(e,)tanhC, (^^Q^ha-cosh^y^.M.E.cosh,. 



\cosh^i -I- cosh ^2 



(33) 



Here x < denotes the distance from the boundary. The full series can be written as Fredholm determinanlj^S. 
Using ([26l) . ([27|) and (I28p . we can calculate the expectation value of the spin operator in the boundary bound state 



(lBk^(:r)|lB) =-ao + i^ / ^ m) tanh^e^^^^-^^^ + . . . 

= -ao + i(To r ^K{0 tanhee2*^"™"'^« + ao tan-, 
Jo 27r 2 



(34) 



.2Mi 



where the contour 7 is defined in Fig.|3]and we recall that k = 1 — h? /2M = cosv. We observe that (1b| (^"^{x) |1b) 
— tTo for X — > — cxD as expected. In the limit w — > 7r/2 the last term remains comparable to cto even for large distances 
X to the boundary. In this limit the boundary bound state becomes weakly bound and its size diverges. The local 
magnetisations of the ground state and the boundary bound state are shown in Fig. [S] For free boundary conditions 
the two results equal each other up to a global minus sign. For finite boundary magnetic field the spins in the vicinity 
of the boundary are aligned parallel to which increases the local magnetisation. 



VI. SPIN-SPIN CORRELATION FUNCTION 



We now turn to the calculation the leading terms in an expansion in powers of the boundary reflection matrix K of 
the two-point function of Ising spins (|17p . Using a spectral representation of the correlator in terms of the scattering 
states ([5]) we obtain the following expression 



00 00 

C(t, Xi,X2) = X! X! ^2 n 2m 



(35) 



where 

C2n2m[T,Xi,X2) ^—rj^ , / , K {^i} . . . K (^^) 

ml (2n)! Jq (27r)'" (27r)^" 

(36) 

(0| a^{T, Xi) 16*1, ... , e2n) {92n, . . . , 6*1 1 (7^ (0, Xa) | -^1 , Cl , ■ • ■ , -Cm, Cm), Xi<X2, 
{0\a^{0,X2)\0i,...,92n) (6'2n,...,6'lk^(T,a;i)|-Cl,Cl,---,-Cm,Cm), a;i > X2. 

We label the various terms in the double expansion (I35p by the numbers of particles in the intermediate state 2n and 
in the boundary state 2m, respectively. The connected correlator is given by 

Ccon„(T,a;i,a:2) = {Ob\% Sa^{T, x^) Sa^^O, X2) \0b) = C{t, xi, X2) ~ {OB\<y'ixi) \0b) (Osk^lxa) |0b) , (37) 
where (5ct^(t, x) = o-^{t, x) — (Ob| cr^(T, x) |0b) and the second term is in fact independent of r. 

A. Regularisation 

The functions p6p contain matrix elements of the form 

(0i,...,0„|a^|Ci,...,Cm), (38) 

which possess kinematical poles whenever 9i — and therefore need to be regularised^. Let A denote a set of one- 
particle excitations and Ai and A2 a partition of A. The scattering matrix arising from the commutations necessary to 
rewrite \A) as IA2A1) is denoted by SaAi, *-e-, |^) — SaAi \A2A1) = Saa2 |^i^2)- For example, if |^) — \6i, . . . , 6*5) 
and \Ai) = {02, O3), then 

\0i,..., 65) ^ SAA^\01,0i.05, 02,03) with Saa.^I- (39) 
If A and B denote two sets of one-particle excitations, the regularisation of the form factors p8|) reads^ 

{A\a'\B) = d{B2)SAA,SB,B {A2\B2) {Ai+iQ\a'\B,) (40) 

A=AiUA2 
B=BiUB2 

= SAA2SB2B {A2\B2) {A,- iQ\ a' \Bi), (41) 

A=AiUA2 
B=BiUB2 

where the sums are over all possible ways to break the sets A = AiU A2 and B — Bi U B2 into subsets. The scalar 
products (A2I B2) are easily evaluated using ([6]). The factor d{A) is present by virtue of the semi- locality of the spin 
operator with respect to the fundamental field and is given by 

d{A) = (-1)"(^), (42) 

where n{A) denotes the number of elements in A. As all rapidities in the remaining matrix elements are distinct, they 
can be evaluated using the crossing relations 

(0,,±iO,...,0,^±iO|a=^|Oi,...,Cj,) = (O|a^|0.,+i7r±i7?,,,...,0,^+i^±i77,^,C,,,...,Cj,) ^^^^ 

We stress that equations (PO]) . (|¥T|l and are valid only for the operator in the Ising model. For the disorder 
operator /i^ we have /^^^^ = 1 and as a result the factor d{B2) in (|40p needs to be dropped, whereas (|¥T|) remains 
unchanged. In other theories and in particular for operators with Iq^i 7^ ±1 (where ^E* denotes the fundamental field), 
additional phase factors related to the non-locality of the operators appear. 

We have checked the validity of (|40p and (|¥T|) using a finite-size regularisation for the Ising model (see below). 
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FIG. 7: Graphical representation of Coo and C2n0- The circles represent the operators , where the upper one corresponds to 
the later time x. The particles created by the A\0)'s are represented by the arrows. We order the arrows between the circles 
from left to right as 6\, . . . , 92n- 



B. Spin-spin correlation function in the bulk system 

The simplest terms in the expansion (|35p are those with m = 0, i.e., terms without boundary contributions. We 
recall the definition of the centre-of-mass coordinates R = {xi + X2)/'2, < and r = X2 ~ xi. We find Coo(t, r) — (Tq (x 
M^/^ as well as 

C2„o(r,r)^-i- r ■ ■ f ^" I /(gi , . . . , r e-^'^"-' ^- -"^^ ^' ^^"^ . (44) 

(2n)! (27r)-^" 

We note that the exponential factor e"^^'''' coshSi ensures the convergence of the integrals. This is a general feature 
of the expressions we obtain. As we will see later, the results for the spectral function remain finite even in the limit 
r — *■ 0. In order to facilitate the comparison to the standard bulk results we now shift the contours of integration in 
as follows. If we write 9 = s + iLp, the exponential factors are given by 

^ — M\r \ cosh gi sgn(r) Mr sinh 9 (sgn(r) Mr cos cp — A/|r | sin ip) sinh s^— (sgn(r) Mr sin i^+M|r | cos (p) cosh s (45) 

which vanishes exponentially for s —t ±00 provided that 

sgn(r) Afr sin<^ + M|r| cosi^ > 0. (46) 

Hence, in the first and third quadrants of the (r, r)-plane, i.e., when rr > 0, the contributions at $He6'i = ±00 vanish 
as long as < (/J < 7r/2 and we may shift the contour of integration 6*,; 9i + i7r/2 without changing the result. In 
the second and fourth quadrants (rr < 0) we may shift 6i ^ 9i ~ it: /2 instead. As the form factors have no poles in 
the strip — 7r/2 < 3m 6',; < 7r/2, we obtain (see also Appendix 

Co nfrr)--^ r TT tanh^^^— ^e-'^^l'^l^-"'"''^'e~*^l^l^''^°'^^^' (47) 

O2„oir,rj - (2^)2" 11 t^nh ^ e e . (4^) 

In order to arrive a this expression we have replaced 9i —9i in the second and fourth quadrants. The result (|47p is 
the well-known result for the spin-spin correlation function in the scaling limit of the Ising model^i. In comparison 
to (|44p we have effectively reversed the Euclidean rotation done in Sec. Illli i.e., we have interchanged t and r back 
in order to interpret r as space and r = i< as Euclidean time. 

If we shift the contours of integration by 6*^ 6*^ ± i arctan (r/r) , ([33]) can be cast in the following form 

i<j 

In real space, t = it, we observe oscillating behaviour for time-like separations and damped behaviour for space-like 
separations. 



C. First-order contributions in the boundary reflection matrix 



In this subsection we calculate the leading contribution in the boundary reflection matrix K. The first, time- 
independent, term is given by 

Co2(t,xi,X2) = -ial r ^K{0 tanhee'*^"^'^"("i'"^)™^h«. (49) 
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FIG. 8: Graphical representation of C02 and C22- The arrows starting at the lower edge represent the particles A^{—^) and j4^(^) 
coming from the boundary state. If these arrows pass the lower circle, which represents an operator a^, it indicates that some 
of the internal particles created by the A^(Si)'s have been contracted with the external lines, i.e., terms like {Oi\ ^) — 2Ti5{6i~ 
appear in the corresponding formulas. The upper index denotes the number of lines connecting the two operators. 



We note that C02 is cancelled in the connected correlation function ([37|) . The first term containing a contribution of 
the boundary to the two-particle continuum is also the first term which contains matrix elements of the form ()38|) . 

2 7n 2tt [2tiY (50) 

g2M max(a;i ,2:2) cosh ^ g — 7\/|r| cosh 9i sgn(r) Mr sinh 9i 

For the evaluation of the second matrix element we can use either ((5D|) or (HI]), which give 

(^^2,elk'h?,0 - (27r)Vo((5(^^l + e)(5(^^2- 0-5(^1 -05(^2 + 0) 

+2^<5((?i - f{02 + iTT + i?72, -e) - 2^(5(01 + f{e2 + iTT + iry2,0 
-2^5{e2 - C) /(^i + i ^ + i ^1 , -e) + 27r5(02 + /(^^i + i ^ + i , 
+.f{d2 + iTT + i?72, 6ii + i TT + i ryi, -e, (51) 

= (27r)2 ao ((5(01 + '5(^^2 - - ^(^i - 5(^2 + 0) 

-27r<5(f?i - C) 7(02 + i^ - i?72, -0 + 2^<5(0i + /(f?2 + i vr - i 772, ^) 
+2^5{e2 - ?) /(^i + i ^ - i ?7i , -?) - 2^5(02 + /(^^i + i 7^ - i ^1 , 
+/(^2 + i^ - 17,2, 01 + iTT - i77i, (52) 

Both regularisations give altogether six contributions, which can be represented graphically as shown in Fig. [HI Either 
regularisation scheme gives rise to a contribution 

C°^ = -ialj^ ^ii'(^) tanh^e2*^™(^i'^=)™"'^«, (53) 

where the additional upper index denotes the number of lines connecting the two operators. The contribution C22 is 
similar to C02 and like the latter cancels in the connected correlation function. 

The next term, C22, actually consists of four contributions which can be combined using the relation K{—^) — 

-m) 

C\2^Tol r ——K{f) tanh ^ COth ^ + ^ ^ g2MficoshC .gn(.) Afrsinhg ^-MKIcoshg sgn(.) Af r .inh ^ (54) 

° 27r 27r 2 2 ^ ' 

For free boundary conditions the integral over ^ has to be understood as a principal value integration. Furthermore, 
the upper sign corresponds to the regularisation ([5T|) and the lower sign to (|5^ . The difference between the two 
is compensated by analogous differences in the contribution C22 discussed below. In order to take the limit 77 ^ 
we shift the contour of the 0-integration in the same way as we did for the bulk terms C2n0- As the exponential 
factors containing 9 in ((54)l equal those in (|44|) . we obtain again the condition (|46)) for the vanishing of the integrals 
at '•KtO — > ±00. Hence, if the space-time coordinates (t, r) lie in the first or third quadrants, rr > 0, we have to shift 
6-^6 -I- i7r/2. In order to avoid the appearance of extra terms arising from the residues of the coth in we choose 
the regularisation ([5T|) . for which the pole in the 0-plane is located at 9 = — ^ — 177, i.e., below the real axis. By the 
same reasoning we choose the regularisation ([52|) if rr < 0. Doing so we obtain 



^1 2 f°° d^dO , i^B^YTxl2 , e + «'±i7r/2±i'7 
^2^2 = Tol J 2^ 2^ ^(e) tanh coth ^ '- 

^^2 MR cosh ^ sgn(r) Mr sinh ^ i sgn(r) Air sinh 6 ^—M\t\ cosh 9 (55) 
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_ _|_ 2 f°° '^^ K(S) Sinh^ 2Afflcoshg -isgn(T) Mr sinh 9 isgn(r)Afrsinh ;^-Af|r| cosh 61 /rg\ 

° 27r 27r ^'^'^ cosh6' T i sinh^ ' ^ ^ 

where we have taken the Umit rj ^ 0. In the second and fourth quadrants of the (r, r)-plane we now change variables 
C ^ and 6* ^ -6, which in particular changes ei''S"('') to I ""I Hence, as i? < we can shift 

^ ^ + i7r/2 and obtain 



CUt,Xi,X2) = -CtI r ^^KU + i^) "^"'^^ gi2Mi^sinh4g-iM|r|sinh9g-M|r|(cosh.+coshO (57) 

+e(/ic-/^) 2ct^ cotw tan^e^^^^^^^e-*'^!^!™" 

rf6> C0Sh6> — Sinu ^^j sinhe ^-A/|T|coshe j-gg-j 

J 27r cosh 6* + sin v 

Here (|58p originates from the pole of K{(^) at ^ = it; (we recall k = 1 — h? /2M = cosw) and is present ii h < he- In 
particular, for free boundary conditions we obtain erg e^*^^ Ko{M\/r'^ + r^) /n. Exphcit expressions for K{£,+i^) 
are given in (jA9p - (|Alip . We will see below that l|57|) yields an oscillating contribution to the spectral function. 

Going through the same steps as above, we find that the final term in the regularisation ([5T|l or ([52|l respectively 
can be cast in the form 

CI. = -4r^r ^^(?)tanhencothi^^cothi±^ 
2 Jo 2TT (2^)2 IJ 2 2 

. tanh^ ^1 " g2M inax(2;i,2;2) cosh 5 g-A.f|r| J], cosh 9. gisgn(r) Mr sinhfli ^gg-j 



r d(hd92 ^2 ^1-^2 A cosh^±isinhf 



2 Jo ^TT J_^ (27r)2 ^ 2 -IJ- cosh^Ti sinh6', 



g2A//max(£Ci,X2) cosh^ g=pi A/|r| X),- sinhfli g— Af |r | cosh Si (60) 

Here we have shifted 6i 9i + iTT/2 for rr > and di 6*^ — i7r/2 for rr < 0, respectively, and taken the limits 
rji 0. Finally, we can substitute 9i — > —di if rr < 0. 
Puting everything together, the result for C22 reads 

C22ir,xi,X2) = -ial /"^if(0 tanhee2^^™'^(^^'^^)™^'^« 
Jo 2^ 

_ 2 /" rf^? j^f^ , . COSh^ — COsh^ „i2Af_R.sinh£ „-i Af k| smh „-A/| r | (cosh 6>+cosh 

y_^2^2^^^^ + '2j cosh^ + cosh^'' 

2 7o 2^7_^ (2^)2 ^ 2 licoshf-isinh^, ^ ^ 

g2A// max(xi,a;2) cosh ^ ^-i Af|r| J2i sinh ^-M\t\ J2i cosh Oi 

+e(/ic - M 2ag cotw tan|e2*^^™"''e~*'^l^l""'' 

d£ cosh ^ - sin V ^_iM|r|sinhe g-A/|r|coshe 

, 2tt cosh 6 + sin v 

The last term is the contribution due to the boundary bound state. 

If one uses regularisation (|52p instead of (|5ip in the first or third quadrants of the (t, r)-plane, one has to keep 
track of the poles of the coth in as well as (|60p in the strip < Urn 6* < 7r/2. A straightforward calculation shows 
that the contributions of these residues cancel out and that one finds again the result (pT|) . In order to check the 
regularisation scheme, we have performed a finite-size regularisation of C'22- The results are presented in Appendix [Bl 
and equal ([6T|) within a relative error of less than 10~^. 

In order to study the light-cone effect in more detail, let us reconsider C22 in the region < h. We start with ([5^ 
and shift the contours of integration 9 ^ 9—i9o and ^ ^ ■C+i'fo; where — arctan(|r|/|r|) and ^0 — arctan(|r|/2|i?|) . 
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FIG. 9: Graphical representation of C24. The first term, C24^, is completely disconnected. In the second term, C24'', the 
contractions of the particles in the intermediate state and the particles in the boundary state are crossed and the operators 
become intertwined. 



This yields 



C',,= all ^^K(^ + i^o) ^°^h(g-i^o)+isinh(e+i^o; 
22 oy-oo27r27r cosh(6'-i6'o) -i sinh(e+iCo) 

We observe the same qualitative features as for the Green's function ([31]) . i.e., oscillating behaviour for < as 
well as AB? < . 

In Appendix [Cl we calculate further contributions to first order in the boundary reflection matrix, but with higher 
numbers of particles in the intermediate state. In particular, we show that C|2(t, xi,X2) — C22{t,X2,xi). 



D. Second-order contributions in the boundary reflection matrix 

In this subsection we calculate the leading contributions to second order in the boundary reflection matrix K. The 
first term of this kind, C04, drops out in the calculation of the connected correlation function. The next term is C24. 
For simplicity, we will restrict ourselves to the first quadrant in the (r, r)-plane, i.e., r > and r = X2 — xi > 0: 



1 f°° d^ld(,2 f°° deid92 ^^^^^ ^ ^^1^ ^2MX2 r.. coshf. f,-MrT., cosh 6, 



C2a{t,x^,X2) = - -7^mi)K{(,2)e 



4 J, (27r)2 (2^)2 



(63) 



(64) 



Using (PO)) the second matrix element can be written as 

(02,0ikN-Ci,ei,-6,6) = (02,0ii-ei,ei)/(-6,6) + (02,0ii-6,e2)/(-ei,ei) 

- (02, 0ii Cl, 6> /(-a, -6) + (02, 0ii - a, 6) /(a, -a) 
+ (02, 0ii a, -a) /(-a, a) - (02, 0ii - a, -a) /(a, a) 

- (02i a) /(0i + + ir,i, -^1, ^1, -^2) + (02i - a) /(0i + + i77i, -a, a, a) 

- (02 1 a) /(0i + iTT + 1771, -a, -a, a) + (02I - a) /(0i + itt + ir^i, a, -a, a) 
+ (0ii a) /(02 + iTT + ir,2, -a, a, -a) - (0ii - a) /(02 + itt + 17^2, -a, a, a) 
+ (01 1 a) /(02 + iTT + 1772, -a, -a, a) - (0ii - a) /(02 + itt + i 772,a, -a,a) 

+ /(02 +i7r + i7;2,0i +i7r + i77i,-^i,^i,-^2,a)- 

A graphical representation for the contribution (|63p using the decomposition (|64p is shown in Fig. [51 We observe that 
there exist two different terms in which the operators are not directly connected. The first one, C24^, is completely 
disconnected and hence does not contribute to the connected correlation function (j37p . The second one, 02^, is 
obtained from the second and third lines of ([M[l . Although there is no line connecting the two operators, the 
contractions of the particles in the intermediate state and the particles in the boundary state are crossed and the 
operators become intertwined. Hence, this term does contribute to the connected correlation function. It is explicitly 
given by (for free boundary conditions the integrals have to be understood as principal value integrations) 

/-fO,b/ ^ ^ \ "'O d£,ld^2 r^,^ s s 2 a^^a 2A/i?5:;, cosh^i ^iMrV;, sinh^, 

C24 (T,a;i,a;2j = — - J -^^-p- A (fi) K (^j tanh — - — e e 
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-e(/ic - h) 2al cot f tan ^ gSMi? cos g-M|r| sin « 

. r ^K{^ + i^) tanh^ ^ + ^ ^""/^ ~ ''^ gi2Mfisinhg ^-A/.cosh;^ (gg) 

We note that the result for free boundary conditions is obtained as the hmit v ^ 0. In the derivation of ([SS]) we used 
that the function had no pole at ^1=^2=11^- Contribution (|65[) can be rewritten as 



-4 r §^^(ei+ieo)i^(6+ieo)tanh^^e-^^^^^^^'^-^^ (67) 

2 J-00 v-T^r 2 

which again shows oscillating behaviour for AE? < in real time. The other terms appearing in (|64p yield sub-leading 
contributions to the correlation function and are determined in Appendix [Cl 

The calculation of higher-order terms in the boundary reflection matrix can be performed along the same lines as 
above (see Appendix [C]) . Although no principal problems appear, the calculations become soon rather tedious. The 
expectation is, however, that the first few orders in K will give accurate results even rather close to the boundary. We 
show below that this is indeed the case for the local spectral function, where the terms calculated above are found to 
be sufhcient for > 0.2. The calculation of terms in the spectral representation with a higher number of particles 

in the intermediate state is also discussed in Appendix [Cl The corresponding terms in the local spectral function are 
found to be negligible compared to the two-particle contributions. The analogous behaviour for the bulk Ising model 
is well documented^i^^'^, and it was argued by Cardy and Mussardo that this behaviour is a general feature of form 
factor expansions in integrable field theories^. 



E. Limiting case 

Let us consider the behaviour of the two-point function in the limit R — > —00 in more detail. We have min(a;i, X2) 
—00, but max(a;i,X2) may remain close to the boundary. Hence, in addition to oscillating terms like C22 there are 
contributions of terms like C22, (^42 ^^^d CI4. These terms diverge in the limit max(xi,X2) — *■ 0, i.e., the series 
expansion (j35p ceases to converge. In the case min(a;i,a;2) = one has to deal with an operator located at the 
boundary. The properties of such boundary operators and their counterparts in the bulk are, in general, very different. 
For example, has the boundary scaling dimension 1/2, while its scaling dimensions in the bulk are (1/16, 1/16)^4 
The treatment of boundary operators in the framework of integrable field theories and form factor expansions has 
been put forward recently by Bajnok, Palla and Takacs^. 

The scaling behaviour of the two-point function in the case of free boundary conditions was derived for the lattice 
Ising model by Bariev"^^. He obtained the following result in the region y/r^ + t2 > 2i? > 1/M, 

{0B\6a^{r,x,)6a\0,X2)\0B) - (^2+^^2)1/4 e''''^ e''^^^^ . (68) 

In the field theoretical calculation the leading term of the left-hand side equals C22, which is according to ([57)1 and 
((58)1 given by 

^12(^,^1,^2) - ^e2^^ifo(Mv/^^^+^) (69) 

,^Tl\ COsh^ — COshff ^i2Af_Rsinh; ^-i Af|r| sinh 6 ^-M|r | (cosh g+cosh g) /jqn 

,27r27r ^ ^ ^ cosh ^ + cosh 61 ' ^ ' 



In the region Vr^ + r^ 3> 1 /M the asymptotic behaviour— of the first term is given by ([68]) . On the other hand we 
can evaluate the second term in the stationary phase approximation, which yields 



"° ■^(go + if) + \ , ^ :^2 , _2.M L , _2UM . Co^iarctanf^). (71) 



2ttM 2/ ^^^^+^4^2+^2 (4^2+^2)1/4 (^2+^2)1/4' 



The condition ^/r^ + > 2R ^ 1/Af translates into > 4x1X2 > 0, where the second inequality follows from the 
fact that the result for C22 is valid only if both operators are located away from the boundary. The condition > 0, 
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however, implies that (|7ip is exponentially suppressed compared to ((69)) . Hence the asymptotic behaviour of C22 is 
given by (|69p and equals the result derived by Bariev in the lattice Ising model. 

We would also like to comment on the conformal limit M ^ of the Ising model. The boundary conditions 



compatible with conformal symmetry are the fixed and free ones, 
derive the two-point function of a^. If one considers the limit X2 — 



3/8 



(OBk^(l/M,l/M) (7^(0, 0:2) |0b) 



Using the method of mirror image s 1 one can 
while keeping xi and r fixed, one finds 

free boundary conditions, 
fixed boundary conditions. 



(72) 



Although we cannot apply the expansion (|35p in this limit, we mention that this result is in agreement with the 
leading correction to the correlation function derived above. For min(a:i,a:2) and t kept fixed this correction is given 
by ([49]) , which is negative in the case of free boundary conditions and hence suggesting the vanishing of the correlation 
function, but positive for fixed boundary conditions, thus supporting a diverging behaviour. However, in order to 
study the vicinity of the conformally invariant point more accurately one has use different methods like the truncated 
conformal space approach2iiS^i2£. 

Finally, we mention that it is not possible to resum higher-order contributions in if by a geometric series, which 
can be achieved for finite-temperature correlation functions in the Ising model^ and the non-linear sigma model^. 



VII. SPECTRAL FUNCTION 

In order to gain further physical insight in the spin correlations it is useful to calculate the corresponding spectral 
function. The positive frequency part is given by 



C(w,Xl,X2) 



dr e'""^ C(t, xi,X2] 



(73) 



uj — > — i aj+(5 



where the analytic continuation of the frequencies is Co — ~iuj + 6. 

There are essentially four different types of terms in the series (|37|) . which we will discuss below separately. First, 
we recover the known results for the bulk from all terms with to = 0. Second, there exist terms which essentially yield 
corrections to the two- and four-particle continua already present in the bulk. The leading corrections of this kind 
are C22 and C42 as well as CI2 and Cq2- Third, we find oscillatory behaviour, which is present even deep in the bulk 
(M|i?| ^ 1). Here the leading contributions are C22 and C24 ■ Fourth, we discuss the contributions of the boundary 
bound state for boundary magnetic fields h < he- The numerical evaluation of the integrals appearing in the spectral 
functions was performed using the VEGAS routine for Monte Carlo integration^. 



A. Bulk Result 



As we have mentioned before, the bulk results for the n-particle level are given by C2no- The full correlation 
function in the bulk is 



00 2 ro 



dOi . . . dd2n 



2n 



tanh 



2 p-iM\r\ Y,i sinhfli -M|t| Y.i cosh 6 



(74) 



i<j 



If we perform the Fourier transformation r — s- a; and calculate the spectral function S'buik = ^nxCbuik, we find for the 
two-particle continuum 



all, d0,de2 , ^2 di 



2 / f">-^^^ 



tanh 



5 ( A/r . sinh 9, ^ s(u; ~ mJ^- cosh 9, 



(27r)2 2 

where we have restricted ourselves to positive frequencies u> > 0. Resolving the (5-functions gives 



(75) 



Arcosh ( 



2MJ 



'5bulk(w,r-) = 2cr^ 



d9 tanh^ 9 

y/tu'^ _ 4Af2 cosh^ 9 



cos (^\r\\/ uj^ ~ 4M2 cosh^ 9 



(76) 
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FIG. 10: Total two-particle contribution to the local spectral density (179 
The bulk result is obtained for MR —oo. We observe that for MR = 
small. 



for fixed boundary conditions (in units of ag/M). 
-0.1 the second order contribution is already very 



We stress that ()76|) vanishes very slowly for large separation r at fixed frequency uj. A similar behaviour is found for 
the two-point function of the disorder operator whose leading contribution is given by the one-particle peak (see 
Sec lVIIl]) . The four-particle continuum can be readily calculated from the corresponding term in (j74p . One finds that 
the ratio S'^uik/'S'buik ~ 1/150 at uj/M = 25 and smaller for lower energies. This suppression of the higher-order 
terms in the series ([74| is a well-known feature in the Ising model and other massive theories^!^'^'^'^'^. 



B. Two- and four-particle continua 



We next turn to terms in the spectral representation (|35p . which have a direct counterpart in the bulk (|74p . 

The two-particle continuum to first order in the boundary reflection matrix is given by C22 and C42, i.e., terms 
with two lines connecting the two operators. Performing the Fourier transform, analytically continuing and taking 
the imaginary part we find after a straightforward calculation using the formulas given in Appendix [Xl 



de 



K{^) tanh^ 



g2A/i?,cosh5 



tanh 



where 



27r J_g, 2tt ^[cj - A/ cosh 61)2 - cosh^ ^ + sinh^ 9 cosh^ ^ + sinh^ 9 

(cosh^ ^ — sinh'^ 9) (cosh^ ^ — sinh^ ^) — 4 cosh^ ^ sinh 9 sinh 9 

■ cos (M|r| (sinh + sinh0)) cosh(A'/|r| cosh^) 
f 2 cosh £, (cosh^ £, - sinh^ 9) sinh 9 + (cosh^ ^ - sinh^ 9) sinh 

• sin(M|r|(sinh6' + sinh^)) sinh(M|r| cosh^) L 



f UJ \ f UJ 

= Arcosh cosh 9 and 9 = Arcosh 1 

\M J \ M 



(77) 



(78) 



We observe that the integrand is exponentially suppressed for large distances from the boundary by the factor 
g2Mflcoshj_ Furthermore, we note that ([771) particle- hole symmetric. The two-particle continuum to second order 
in K, which we denote by 5*1^, is obtained from Cl^, C^'l and CI4. Its explicit form is derived in Appendix ID II 

The total two-particle contribution to the local spectral density is given by summation of all terms with two lines 
connecting the two operators. In Fig. [TO] we have plotted the truncation at the second order in K, 



SI,{lo, R,r^O)^ 5^,ik(co, r = 0) + 5iV(^, R,r = 0) + SIk{uj, R,r = 0) + 



(79) 



for fixed boundary conditions. The results show that the effects of the boundary can be treated in perturbation theory 
for M\R\ > 0.2. This is also supported by the study of the ratio of to Sfj^ at their respective peaks, which is 
plotted in Fig.[TT| For example, for fixed boundary conditions and MR = —0.1 the latter is about 5%. 
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FIG. 11: Ratio of the first and second order (in the boundary reflection matrix) contributions to the two-particle continua 
S2K/S1K fl't their respective peaks for different distances from the boundary. 



Starting from C42 and Cq2 we have calculated the four-particle contribution to the spectral function Sfj^ to first order 
in K. The explicit result is given in Appendix ID II Comparison with Sfj^ shows that the four-particle contribution 
is negligible in the low-energy region lu < 25M. This indicates that the fast rate of convergence in the number of 
particles in the intermediate state, which is a well-known fact for the Ising model in the bulkl i^^i'^^i'^^ , is also present 
in the system with boundary. 



C. Oscillating contributions 

The corrections to the bulk we have studied so far are negligible for AIR <ti —1. In this section we discuss oscillating 
terms which affect the spectral function also deep in the bulk. The leading term of this kind is given by (|57p , we will 
denote it by 022^"^ ■ The corresponding contribution to the spectral function reads 

Arcosli ( -fj — 1 ) 

:,l,osc/ , n „ _ n^ _ ^'^Q 



C.i,osc, r, n^ ^Cr^ f d( W - 2M COsh ^ 

S22 R,r^O) 



27r ^(w-Mcosh^)2-M2 (gQ) 
meK{^ + i^) cos(2Afi?sinh^) -3mX(^ + if) sin(2Mi?sinh^)j , 

where K(^^ + i f ) given by (|A9[) - (|Alip . For simplicity we have restricted ourselves to r = 0. We stress that ([M]) is 
not exponentially suppressed for large distances from the boundary. In fact, the integral ([50)1 possesses a square-root 
singularity at the upper limit ^ = Arcosh (-p^ — l). Expanding the non-singular part of the integrand about this value 
yields 

ci.°sc, p n^ co8{AR + ip) j— — — - 

(w,i?,r = 0) — . , A = 2Vt^2 „ 2A/W, (81) 

where depends on the boundary conditions. 

The next term we want to analyse is 0'^^'^"°^'^ . The corresponding spectral function is found to be 



.2 /•?' 



rtO,b.osc/ D \ / ""^ 4 T,: 

^24 {uj,R,r) ^~ — — tanh 



^ 2M J_^, 27T ^(uj - Af coshO^ - 2 

5Hei^ + i f ) ^zK{^ + i f ) - ^mK{^ + i f ) 3mK{i + i f ) 1 cos(2A/i?(sinh^ + sinh^)) (82) 

- 2 D\eK{^ + i f ) 3mX(C + i f ) sin(2Mi?(sinhC + sinh^) 1, 
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FIG. 12: Oscillating terms 522°"'^ (black line) and 824''°^'^ (red line) as well as S^'^"^'^ + S^^f^'^ (green line) for fixed boundary 
conditions, uj = 5M and r = (in units of oq/M). We stress that £'22°'^'^ and 524''°°'^ are present deep in the bulk, whereas 
524°'^'^ + ^44°'''^ is strongly suppressed for large M|J?|. 



where ^ and ^' are given by 

I = = Arcosh - coshe) , i' = Arcosh - l) • (83) 

We note that ([5^ is independent of r. Numerical evaluation of the integral in the range MR < —20 and ui « 5M 
shows that 

oO.b.osc, p „s C0S(A^.R + if') . , . 

524' (c^,i?,r = 0) ^M\j^\y/2 ' ^>^- (84) 

The larger oscillation frequency A' is due to the additional term sinh^ > in the sine and cosine in ([5^ . Although 
(151)1 falls off faster than (|8ip . 5*24'^'°'''^ should not be interpreted as a higher-order correction to 522°*"^, as both terms 
are present deep in the bulk. The analysis of the corresponding contributions in the presence of a boundary bound 
state (see below) rather suggests that 5*24^'°'"^ represents an independent oscillating contribution of order K; and hence 
that the total oscillating contribution to first order in the boundary reflection matrix is given by 

Sr^icj, R, r) = 52V(^, r) + Sl'^^°'''{u,, i?, r). (85) 

The higher-order corrections in X to ((501) and ^ are in fact obtained starting from Cl'^'"'' + C];^'-'' and Cl^""" + Cl^""" , 
respectively. The corresponding spectral function 5*2!^ is evaluated in Appendix ID 21 We note that the appearing 
integrands are exponentially suppressed for large distances from the boundary. 

We have plotted ([80)l and ([82)1 as well as (|D5[) in Fig. [12] for fixed boundary conditions, fixed energy and varying 
distance from the boundary. The first two are present in the bulk [MR ^ —10), whereas the third one is strongly 
suppressed for large M|i?|. This underlines the interpretation of (jPSp as first correction to ([50)1 . Finally, the correction 
to ([50)1 with three particles in the intermediate state is given by 6*42°'^'^. It is found to be negligible for all values of 
MR and energies u) < 25M. The same result holds for C|4 the first correction to (|82p with more internal lines. 
This again underlines that the series ([55)1 is rapidly convergent as n ^ 00. 



D. Contributions from the boundary bound state 



In the previous section we have investigated the oscillatory behaviour of the two-point function. This behaviour 
was caused by the terms C22, C24'', C24 -I- C^^ and + C^g. In Sec. IVII we have seen that in the range h < he all 
these terms possess additional contributions due to the existence of the boundary bound state, which we will now 
discuss in detail. 

We first consider (j58p . which we will denote by C2^^^. The corresponding contribution to the spectral function can 
be cast in the form 

ii bbs/ n 2crQCotwtan| uj — 2Msmv 



Slt%u;,R,r) = ,/ " ; :;^ = cos(rV(^-Afsin.)2-M2) e^*^«™«^ (86) 

y/{uj — M smv)^ — M^ ^ 
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and is seen to be exponentially suppressed in MR. We recall that the parameter v is defined as k = 1 — h? /2M = cos z;. 
In the limit w — > 7r/2 the suppression of 5*22^'^*' becomes less and less effective. In this limit the boundary bound state 
becomes weakly bound and its effective size diverges^i. 
Similarly we obtain 



/-rO,b,bbs/ 7-) \ 
624' (t^,it,r) 



2KPal cotv tan^i ^2MRcosv 



M svav) 



A/2 



-5nei^(C + if) (cosh2|-cos2w- 2) -23mX(^" + if) sinh|cost; 



cos(2Afi?sinh^~) 



(87) 



sinh^ cost! H — 3mK{i + i^) (cosh 2^" — cos 2v — 2] 



sin(2Afi?sinh|) 



where 



^ — ^{uj,v) — Arcosh y — — sinv 



(88) 



The contribution 824^'^^^ is independent of r and oscillates in R with the frequency-dependent wave number 
2-\/ (cj — Af sinti)2 — M^. In particular, for large ui this wave number approaches 2{oj — Msinv). Furthermore, 
(j87p is exponentially suppressed for i? ^ — 00. 

Interestingly, although both (|86|) and ((87)) have threshold singularities at w = Af + Af sinw, these singularities cancel 
each other. This shows that we have to interpret their sum as the first-order contribution of the boundary bound 
state to the spectral function, i.e., we have 



S^'"K^{u;,R,r) = Sli'''%u;,R,r) 



Si 



■O.b.bbs 



iLO,R,r). 



(89) 



This indicates as well that the similar sum ([SS]) should be viewed as the oscillating contribution to first order in K. 
We observe that has a gap M-l- Afsinu which equals the energy of the boundary bound state plus one additional 
particle. On the other hand the spectral function vanishes at w = Af sinw. Combining these two observations we 
conclude that 



(lBk='|OB)=0 but (lB,0k^|OB) 7^0, 



(90) 



where we have used the notations in the original picture with the boundary located in space. The vanishing of the 
spectral function for w — > Af -t- Af sin v means that the rapidity-dependence of the matrix elements of overcomes 
the van-Hove type singularities in the one-particle density of states at 9 = 0. 

The result for free boundary conditions is obtained in the limit v ^ 0. In this case simplifies to 



S'^^^%,iuj,R,r) 



Jl „2MR 

7q e 



A//2 



cos(rA/a;2 _ ^/2^ _ ^ [ cos{2R^/uj^ - NP) + 



Vlj2 - A/2 



sin(2/? Vtj2 _ 



and the spectral gap equals M. 

The second-order contribution in K of the boundary bound state to the spectral function is given by 



(91) 



52^]f(u;,/?,r) = 5. 



^'^'^''(c., R, r) + 5l^'=^(c., /?, r) -f Sl^^io, R, r) + S^t^u;, R, r). 
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(92) 



Explicit expressions for the various contributions are given in Appendix ID 31 We find that all integrands appearing in 
82^^ are exponentially small for large A/|/?|. Furthermore, S'jjf remains finite at the lower threshold uj M + M sinv, 
which is ensured by the cancellation of all singularities in ([92)1 . We expect similar cancellations of singularities to 
occur at all orders. It is also straightforward to calculate the first correction to ([89|l with more internal lines, which 



1 ,^3,bbs 

IS given by C42 



^2^b,bbs^ Wc find this to be negligible for all values of R and energies oj < 25 A/. 



E. Results for the local spectral function 



In the preceeding sections we have discussed individual terms in the expansion of the spectral function. We are now 
in a position to combine the various contributions and present results for the local spectral function of Ising spins up 
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FIG. 13: Local spectral function (|93|) for MR = —0.5 and fixed and free boundary conditions as well as a boundary magnetic 
field corresponding to u = 7r/6 (in units of a^/M). The spectral function possesses a gap of 2M for fixed boundary conditions. 
For boundary magnetic fields h < he the gap is smaller than 2A4" and goes to M in the limit of free boundary conditions. 




2 4 6 8 10 12 14 16 18 20 22 24 



FIG. 14: Local spectral function (|93p for a boundary magnetic field corresponding to v = n/6 and different distances from the 
boundary (in units of a^/M). 



to second order in the boundary ii'-matrix. If we truncate the expansion at the two-particle level as well, the local 
spectral function truncated is given by 



mcated is given by 

S'^'icj, R) - Sii,^{u;, R) + S\%{lj, R) + Sl%{u;, R) + . . . 
S'^U{^,R) = SL,^{u;,R,r = 0) + ... 



(93) 

S'^Ui^,R) = Sl,,^iu;,R,r = 0) + ... (94) 
S'}f(w,i?) = Sfi^iuj,R,r^O) + S°'j^{uj,R,r = 0) + S^]^%uj,R,r = 0) + ... (95) 
S^i^{uj,R) = SiK{uj,R,r = 0) + S°'i^{uj,R,r^Q) + S^^%Lu,R,r = 0) + ... (96) 

where the terms originating in the boundary bound state are only present if h < h^. The dots in (|93p represent higher 
orders in K, whereas the dots in ||M|) - ([M)) stand for terms with more than two particles in the intermediate state. 
As we have shown above, the truncation at second order in K is sufficient for MR < —0.2. We have further argued 
that the truncation at the two-particle level gives a very good accuracy for energies tu < 25M. 

We have plotted (f93|) in Fig. [13] for fixed and free boundary conditions as well as a boundary magnetic field 
corresponding to k — cos(7r/6). For fixed boundary conditions the local spectral function has a gap 2M. The same 
behaviour of is well-known from the system in the bulk. For sufficient small values of the boundary magnetic field, 
however, the gap is given by Ei + M, where Ei is the energy of the boundary bound state. This shows that the 
matrix element of cr^ between the vacuum state in the presence of the boundary and states containing the boundary 
bound state and one additional particle is non-zero. In the limit of free boundary conditions Ei tends to zero and the 
gap approaches the single-particle gap M. Furthermore, we observe oscillating behaviour in lu for all three boundary 
conditions. 

In Fig. [T3] we plot ([M)) for a boundary magnetic field corresponding to k = cos(7r/6). Close to the boundary the 
spectral weight is concentrated at low energies and in particular at energies u! < 2M. The weight below 2M is due 
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FIG. 15: Local spectral function 
al/M). 



93)) for fixed boundary conditions as function of the distance from the boundary (in units of 



to the creation of the boundary bound state and one additional particle. The term responsible for this behaviour is 
((89)) . Both ()86p and ((87)) are exponentially suppressed with increasing distance from the boundary. Hence the spectral 
weight in the region < w < 2M is very small, as can be seen in Fig. [TDfor MR = — 10. This strong suppression 
may complicate an experimental detection of this effect of the boundary bound state on the local spectral function. 
Furthermore, the spectral function oscillates in oj with a frequency depending on the distance from the boundary. For 
large \R\ these oscillations become very rapid, which may cause difficulties in detecting this effect as well. 

Finally, we plot the local spectral function as a function of the distance from the boundary in Fig. [15] It is clear 
from this plot that the oscillations are only algebraically decaying in \R\. 

In summary, the local spectral function ()93p exhibits two new features in the presence of a boundary. First, the 
spectral function shows oscillatory behaviour both as a function of uj with a frequency depending on the distance 
from the boundary, and as a function of R with an energy dependent frequency. Second, for values of the boundary 
magnetic field below the critical value ~ \/2M we have found spectral weight in the interval AI < uj < 2M, i.e., 
within the gap of the bulk Ising model. This spectral weight is due to the existence of a boundary bound state, and 
hence strongly suppressed with increasing distance from the boundary. 

Any quasi one-dimensional material described by the quantum Ising model needs to be thought of as an ensemble 
of finite length chains. If the average chain length is large, a model in terms of semi-infinite chains should constitute 
a good starting point and the results obtained here should be applicable. In particular, it should be possible to 
detect the midgap states in inelastic neutron scattering experiments. On the other hand, the oscillatory behaviour of 
dynamical spin correlations is likely to average out. 



VIII. CORRELATION FUNCTION OF THE DISORDER FIELD 



In this section we briefly discuss the two-point function of the disorder field /i^. Up to first order in K we find 
(ObI % M^(r, xi) /i^(0, X2) |0b> - (0| % 2:1) ^'(0, X2) |B) 

= (0| fi^ir, xi) /i^(0, X2) |0) + r ^ K{0 (0| m^(t, xi) /i^(0, X2) hC, + ■ ■ • , (97) 

where we have assumed that xi < X2. The first term equals the result in the bulk. The single-particle contribution 
to it is 



de 



de 



J -00 27r 27r tt 

where we have used ()A7I) . After analytic continuation, t = ii, we observe typical light-cone behaviour. Fourier 
transforming we find that the single-particle contribution to the bulk part of the spectral function is given by 



COii{r\/~UJ 



_ M2 

We note that the spectral function shows undamped oscillations as function of r. 



(99) 
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If we assume r > 0, the leading contribution due to the presence of the boundary can be calculated using ([40]) 



(^lM'he,0 = {e±iO\fi^\~t,0+'27r<JoS{9~0~2naoS{9 + 0. (100) 

The connected piece merely yields a correction to the one-particle spectral function (|99p which vanishes for large 
distances from the boundary. However, the disconnected pieces give 

r ^ i^(^) e2Mfl cosh e^i Mr sinh e ^ ^2 r ^ ;^(^ + J cosh « (^q^) 

J-oo 27r 27r 

+Q{h^-h) 9(^0 -w) 2(7^ cott- tan^e^^^^^'^^e-^^^™", (102) 

where = arctan(T/2|i?|). Note that (jlOip and p02p are independent of r. If we perform the analytic continuation 
T = i< in (jlOip . we find an exponentially damped behaviour for 4i?^ > , and an oscillating behaviour for AR^ < . 
The contribution to the spectral function corresponding to (llOip is given by 



,osc 

'-'iK — 



9^eX(Arcoshf-+if) cos(2v/w2 _ ^,/2 i?)-3mi^(Arcosh|^+i f ) sin(2v/w2 - NP R) . (103) 



We observe that (|103p is not damped for large distances \R\. This is due to the dissipationless propagation of a single 
particle to the boundary and back, and a likewise dissipationless reflection off the boundary. On the other hand, l|102p 
represents the contribution of the boundary bound state. It is seen to give rise to a i5-peak contribution to the spectral 
function at the energy lo = Msinv, which implies that (1b| /i^ |0b) ^ 0. Higher order contributions to the two-point 
function of disorder operators can be calculated by the same method we employed above for the spin correlations. 



IX. CONCLUSIONS 



Our main result is the calculation of the dynamical spin-spin correlations in the Ising field theory with a boundary. 
We have derived an expansion in powers of the boundary reflection matrix K and shown that it converges rapidly even 
in the case where both operators in the two-point function are fairly close to the boundary. We have also demonstrated 
that like in the bulk case higher-order terms in the expansion in the number of particles in the intermediate state 
are negligible at low energies. The most notable efi^ect of the boundary is that at sufficiently late times the spin-spin 
correlations show oscillatory behaviour arbitrarily far away from the boundary. As is well known, for small values of 
the boundary magnetic field a boundary bound state exists. This bound state leads to a contribution to the spectral 
function within the gap of the bulk Ising chain. Similar features are also found for the Green's functions of the 
Majorana fermions and the two-point function of the disorder field. 

We have seen that the expansion in powers of the boundary if-matrix breaks down close to the boundary. In 
order to access this regime, other methods are necessary. One possible approach is the truncated conformal space 
approach, which has already been applied successfully to the analogous problem for one-point functions2S. In would 
be interesting to generalise these results to the case of two-point functions and in this way obtain accurate expressions 
for the two-point function for all values of MR. 

The results obtained in the present work have applications not only to the quantum Ising model itself. It is 
well-known that both anisotropic spin-1 Heisenberg chains^^ and the weak-coupling two- leg spin- 1/2 ladder— can be 
described in terms of three and four Ising models respectively. Correlation functions of the staggered components of 
the spin operators in these models are represented by products of two-point functions of the spin and disorder fields 
of the quantum Ising chain. Using the results obtained in the present work one may calculate dynamical correlation 
functions for these systems. 

Finally, having established that the expansion in powers of the boundary iiT-matrix converges even relatively close 
to the boundary, one may apply the method used here to other systems. In a forthcoming work we apply the method 
used in the present work to the calculation of the local tunneling density of states in a one-dimensional charge density 
wave stated. 
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APPENDIX A: USEFUL FORMULAS 



We summarise some formulas, which were frequently used in the derivation of the 6*2^ 2m 's. 

sinh(x ± i7r/2) = ±i cosha;, 

cosh(x ± i7r/2) — ±i sinhx, 

, X — y X + y cosh x — cosh y 

tanh tanh = — , 

2 2 cosh X + cosh y 

. x — y^'nr/2 , x + ?/±i7r/2 cosh a; =F i sinh w 

tanh tanh — 



tanh --y^'^-/^ coth - + y^'-/^ 



cosh X ± i sinh y 

cosh y ± i sinh x 
cosh 2/ =F i sinh x ' 



x — y^n:/2 x + y±nr/2 cosh a; ± i smh y 

coth coth = — . 

2 2 cosh X =F 1 smh y 



(Al) 
(A2) 

(A3) 
(A4) 
(A5) 
(A6) 



li x,y £ M., y > 0, f{9) is analytic in — arctan(a;/?/) < 3m9 < arctan(a;/j/) and at most exponentially growing for 
6 ±oo, then 



i Mx sinhS ^— My cosh 9 



(A7) 



where 6*0 = arctan(a;/y) . In particular, if a e R and f{d) = e"^ then. 



2^ 



^ct6 Mx sinh 6 ^~My cosh 6 



i y + X 



a/2 



27r 



(A8) 



If / possesses poles in the strip — arctan(a;/?/) < 3m6 < a,Tctaii{x / y) additional contributions appear on the right-hand 
side of (ITZl) . 

The explicit form of the boundary reflection matrix K at ^ + i7r/2 is 



i^(e+if) 



l + tanh^^/2 -|_ sinh^ ^ 



tanh^ 2 ^ ) ('^^ ^ ^"^^^ ^) ~ ^^'^^ f ^^^'^^ 



-2i ( tanh | (k^ - sinh^ ^) + Ksinh^ f tanh^ | - 1 



(A9) 



For the special cases of free and fixed boundary conditions this simplifies to 

1 



iffroo(C + if) 
^fixcd(C + if) 



1 + coth^ ^/2 
1 

1 + tanh^ ^/2 



9 £. $ 
1 - coth^ - - 2i coth - 



1 — tanh^ ^ — 2i tanh ^ 



(AlO) 
(All) 



We note that the real part is even under ^ — > — ^ whereas the imaginary part is odd. 

In the calculation of the spectral functions we use that if cj G M, f{9i, O2) = —f{—Oi, —6*2) and / sufSciently well 
behaved, then 



/(0i,< 



w-X)iC0sh6li 
j deid02fi0i,O2)s(Lo-'^cosh0,^ = 0. 



(A12) 
(A13) 
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Furthermore, we use 



-pr cosh ^ + i sinh _ n?(cosh^ ^ - sinh^ Oj) - 4cosh^ g J]^ sinhgj 
^-licoshf-isinh6', ^ Yllicosh^ i + smh^ e,) 



1=1 

2 



-^J-coshf-isinh6'i n?(cosh2 ^ + sinh^ 6*,) 



1=1 



and similar formulas for the higher-order terms. Note that (jA14[) is even under simultaneous reflection 6i ~9i, 
while (jAlsp is odd. 

APPENDIX B: FINITE-SIZE REGULARISATION OF C22 

In order to check the infinite volume regularisation scheme for form factors involving two multiparticle states, we 
evaluate C'22 in the flnite system < t < L. A similar analysis was performed recently for the correlation functions 
in the Ising model at finite temperatures^^. The form factors for in the finite system are^ 

■m n 
i=l j=l 

where the rapidities in the Neveu-Schwarz (NS) and Ramond (R) sector are given by 

NS : AfLsinh^fe = 27rfc, fc e Z + i, (B2) 
R: MLsmh9n^27m, neZ. (B3) 

The function Fm n in (|B1[) is the infinite-volume form factor 



m n 



F„,„(Ci,...,C™|0i,...,0„)=i^(™+")/^Jc7o n tanh^^ n tanh^^niT'^^thSl-^, (B4) 

— l *J = 1 *=lj = l 

where the symbol [ J denotes the floor function, i.e., \_x\ is the largest integer I < x, and the constant as well as the 
leg factors are 

5(L) = l-HO(e-^), ff((?)=3(^)-^==L_ + 0(e-^). (B5) 

VML coshe' 

We note that due to the fact that connects the NS and R sectors of the Hilbert space, the rapidities and 9j 
cannot coincide and therefore no singularities occur in the flnite volume. 

In terms of these form factors the finite volume regularisation of C22 is given by (we assume xi < X2 for simplicity) 

C22 ^ Y1 ^(Cfc)Ns(0k^(T,xi)|7ii,n2)R, R(n2,«ik'(0,a;2)|-fc,fc)Ns (B6) 

feeNS ni,n2eR 
fc>0 

2 

^ fcGNS ni.n2eR 

fc>0 

■~g\^,) tanha tanh^ ^IL^ f[g\0.^^) coth coth ^^i±^ . (B7) 

i=l 

We have evaluated (|B7p and ((6T|) numerically for several values of xi, X2 and 0.2 < r < 1 as well as different values of 
the boundary magnetic field. The two expressions coincide within a relative error of less than 10^**. For the evaluation 
we used L = 50 and N = 350 as cut-off for the momenta, i.e., < k < N and —N < rii < N . The result does not 
change if we use Ramond states as the outer states in (jB6p instead of Neveu-Schwarz ones. 
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FIG. 16: Graphical representation of C42 and C52. The dots represent diagrams with the same number of contractions. 



APPENDIX C: HIGHER-ORDER CORRECTIONS TO ([35 



In this appendix we calculate several higher-order terms in the expansion ([3d)) , which can be compared with the 
leading contributions calculated in Sec. IVII In order to simplify the notations, we will restrict ourselves to t > and 
r > 0. The correlation function for general r and r can be obtained as in Sec. IVII 



Calculation of C42 



The next term in the series (|35|) is given by 



deide2d93d94 



2M max(2: 1,3:2) cosh ^ g— A^lr] 52 . cosh 8; 



(CI) 



matrix element 







•e' 


3gn(r 


) Mr sinh 


«7(0 


1,^2,^3, ^4) (^4, 93,62, 




ourselves to r 


> 


and r > 0, 


we 


can apply 


the regularisation scheme 


(|40)) to evaluate 


(04, 03, 02, 01 


a" 


he,e>-(^4. 


03 1 




+ Itt 


+ i?72, 01 + iTT + ir/i) 




- (^4, 


^21 


-C, 0/(^3 + 


i TT 


-|-i?73,0i + 


i7r + 






+ (^4, 


^il 


-C, 0/(^3 + 


i TT 


+ i?73 , 02 + 


iiT + 






+ (^3, 


^2| 


-e,e)/(^4 + 


i TT 


+ i??4,6'i + 


i7r + 


i?7i) 






^ll 


-e,e)/(^4 + 


i TT 


+ i?74, 02 + 


i7r + 


i'72) 




+ (^2, 


^ll 


-C, 0/(^4 + 


i TT 


+ i??4,6'3 + 


Itt + 







- (6'4|0/(f3 +i7r + i?73,6'2 +i7r + i772,6'i 

+ (6'3| /(04 + ivr + irn, 02 + Itt + i7;2, fi + i^r + ijyi, -f) (C2) 

- (02|O/(04 + i7r + i?74,03 + i7r + i?73,6'i 
+ /(04 + i7r + i?74 , 03 + i7r + i7]3,6'2 
+ {04-0 /(^'s +i7r + iry3,6'2 +i7r + i?72, 

- (6*31 -0/(^4 + ITT + iry4, 6*2 + i7r + i7?2, 

+ (02| - /(04 + iTT + i?74, 03 + ITT + i 773, 

- (0i| -0/(^4 + i7r + i?74, 03 + ITT + i 773, 

+ /(04 + iTT + i?74,03 + iTT + iry3,02 + ITT + ir/2,01 + ivr + i?7i, -^,^), 

This yields three terms, which are graphically represented in Fig. [111 The first one is found to be C|2(t, xi,X2) 
C22{t,X2,xi) for general Xi and X2- Furthermore, the contributions due to lines seven through 14 give 



+ i7r + i77i,-^) 
+ i7r + i77i,-f) 
+ 177 + ir/1, -^) 
+ Itt + ir/2, — 
9i + iTT + ir/i,C) 
9i +i7r + ir/i,e) 

91 +i7r + ir/i,e) 

92 + ITT + ir/2,0 



C12{t,Xi,X2) 



d£_ d0id02d03 
2^ (27r)3 

■ f[ tanh^ ^ 

i,i=l 



K{^) f[ tanh ^— ^ coth ^ 



i=l 



^ g2Mi?cosh{ ^-MrJ2iCosh0i Mr (smh J+X]; smh 6 



(C3) 



25 



6 27r (27r)3 + ' 2 j JU ^osh ^ + cosh 0. 



^ ^ 3 2 (27r)3 cosh6l, + sinu 

3 



Here (IC4[) shows the typical oscihating behaviour and (jCSp is present if h < he 
Finally, the last term of (|C2[) yields 



C42(r,a;i,a;2) =-^^J^ ^ J J2n}^ ^^^^ H ~ 



»<i (C6) 

-j-j- cosh g + i Sinh 0^ ^2A/K2Cosh^ p-iMr^^sinhe, ^-Mr^^coshSi 



cosh^ — i sinh^i 

The next term in the expansion psp possesses a similar leading term, formally we find Cg2(T, xi,X2) — C^2i'''^ ^2, xi 



2. Calculation of C24 

The fourth to seventh line of ([M]) yield 

cosh ^1 — i sinh ^2 cosh ^1 + i sinh cosh ^2 — cosh 
cosh^i + i sinh 1^2 cosh^i — i sinhS cosh ^2 + cosh0 

^g2A./a:2 cosh^i 2Mi?, sinh ^2 g-i A/r sinh 8 g-A/T(cosh 0+cosh ^2) (C7) 



df f°° d0 , , cosh f — cos V cosh — sin w 

^ ' ■i^(^)tanh^- ^ 



2ti y_oo 2vr cosh ^ + cos w cosh + sin w 



cosh ^ + i sinh 
cosh ^ — i sinh 



^2Mx2 sinhf ^-i Afr sinh 6 g-A/r cosh S 



(C8) 



We stress that the integrand in (|C7p is exponentially suppressed for X2 — > —00. Furthermore, the last term in 
yields 

d^idi2 f°° d0id02 , ,^ , 2 ^'i - ^'2 /cosh^i -cosh^^ ^ 



^2 . N ^0 r d£_id^2 r d0id02 -TT r^^.s 



(27r)2 (27r)2 11 2 Vcosh^+cosh^ 

TT COSh^j +i sinh 0j ^2Af:r2 coshfj g-iAfry:,sinhei ^-MtY,, coshOj 

^ ^ cosh ti — i sinh . 

1,1=1 ' ■> 
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FIG. 17; Graphical representation ol C44. 
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FIG. 18: Graphical representation of CI4, C26 and Cjg. 



3. Calculation of C44 



For the evaluation of C44 we have 

1 dCid^2 deide2d03d0i 



(27r)4 



2Mx2 I], cosh^i 



(CIO) 



. g-Afr E, cosh 9, gi Afr -nh 9, f(^0^^0^^0^^ Q^) (^Qi, 03,02, 0l\cJ^ 1-6,6,-6, 6) 



If we regularise the second matrix element according to (HO]) we first obtain (see Fig. [T7|) the completely disconnected 
term C^^(t, xi,X2) = Co4(t, 0:2, a^i). Along the same lines as in the calculation of we find after some algebra 
that C\^(t, xi,X2) = C\^(t, X2,xi). Furthermore, there are two different terms with two lines connecting the two 
operators, first 



C,Ar,x„X2)=-- 1 



tanh 



2 2 
d0id02 -TT ,^ , , ^ -TT cosh 6 — i sinh^i cosh 6 + i sinh^^ 

11^(40 tanh 6 11 cosh6+isinhe, cosh 6 - i sinh 

2 ^1 — 02 2A'/ ail cosh ^1 2A/x2Cosh42 „-i A/r E; smh 9; - Afr E; cosh 9; 



(CU) 



which is not symmetric under 6 6, and second. 



cI:^(t, xi,x2) 



^ ^6^6 ti^i'^^2 

CO (2^)2 (2^)2 



i^(6)i^(6) tanh' 



6-6 



tanh 



2 



T-r cosh 0j + i sinh 
^ ^ cosh 0, — i sinh 6 



2 2 

2Af/?Ei cosh^i ^-i Mr J2i sinh Bi A/t Ei sinh g-AfrEiCosh6 



(C12) 



Here we can again shift 6,2 ^ 6,2 + i7r/2. The result will, however, only yield a negligible correction to the spectral 
function. 



4. Calculation of Ce4, C26 and C46 



The last terms of the series psp we wish to calculate are the leading contributions of Cq4, C26 and C46 sketched in 
Fig. [m The first one is found to be C|4(t, xi, 0:2) = C|4(t, 0:2, xi). 
For the evaluation of the second one we start with 



C2,ir,x„X2) = -J^ 



^2Afa;2 JZi cosh?; 



(27r)3 (27r)2 11" 

Air cosh gi A/r E. sinh 9, ^(^^ ^ ) ^ | | ^ ^ ^ ^ ^ 



(C13) 
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If we regularise the second matrix element according to (|40p and keep only those terms in which the two intermediate 
particles are contracted with two particles possessing different rapidities from the boundary state, e.g. the term 
proportional to (02,^i| ^2,^3): we find 



C°6(r,xi,a;2) - i^ ^ if(a) tanha e 



2Mx2 cosh ^1 



^(^2 + 1 2 ) ^(^3 + 1 2 ) tanh 

^ -pj- COshgi -i Sinha ^,2MRT.l=.,s\nhi, g-AfrELsCoshS, 

\^ cosh a + i sinh ^ 



(C14) 



+e(/ic - /i) 21^2 cotu tan-e2*^^'=°''''e-*^^"'"'' 

|l r |a;f(^,)x(e2 + i|) tanhCie^^^---^^- cosha-cos. 
27r ' '^■^ 



• tanh' 



27r 



cosh a + cos V 



6 + i (7r/2 - i;) cosh a - i sinh^ gi2Afi?si„h6 g-A/rcoshc^ 
2 cosh a + i sinh a 



We label (ICTill by Cl^'' and (ICT5ll by Ca^g In the same way one obtains C^q{t, xi,X2) = C2q{t,X2,xi) 



(C15) 



APPENDIX D: HIGHER-ORDER CORRECTIONS TO THE SPECTRAL FUNCTION 

1. Two- and four-particle continuum 

The second-order contribution to the two-particle continuum is given by CI4 -I- C'^^ + CI4, which after Fourier 
transformation (|73p and analytic continuation reads 



C2\(c., ii, r = 0) - (CI4 + Clt -t- Cl^) (a;, i?, r - 0) 



„2MRY,i cosh^i 



cosh 6*4 -hi 5 



• tanh 



2 Cl — P2 



/ cosh a — cosh a 
\ cosh a + cosh a 



2 2 

n 



cosh a + i sinh 6j 
cosh a — i sinh 0j 



(Dl) 



cosh a — i sinh 9i cosh ^ + i sinh 9i 
\^ cosh a + i sinh 0^ cosh ^ ^ i sinh 9i 

When calculating the spectral function, the contributions proportional to the principal value vanish as the imaginary 
part of the terms in the squared brackets is antisymmetric under 6i — > —9i. After some straightforward algebra we 
obtain 



d^id^2 



-fl- 



K(a) tanh a 



tanh 



2 e-9 



,0 {27r)2 J_g, 2tt 1^1 (cosh^ ^ + sinh^ 9) (cosh^ ^ + sinh^ 9) y/{u - M coshff)^ - 

g2Af _R X); cosh C 2 



(cosh^ a + cosh^ a) TTl(cosh^a -sinh^6')(cosh^a -sinh^e*) -4cosh^6 sinh 6* sinh ( 
(cosh a + cosh 42) I f^i 

-hSjJcosha (cosh^a - sinh^ 6*) sinhe* + (cosh^a - sinh^ 6*) sinh6' i. 



i=l 



(D2) 
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where and 9' are defined in (|78|) . We stress that the integrand is exponentially suppressed for large M\R\ by the 
factor e^^^^LiCosh?.^ 

The first order correction in the boundary reflection matrix to the four-particle continuum is given by + (^62 ■ 
The corresponding spectral function is 

Stt^i^, R,r^Q)= + St^) {io, R,r = 0) 

al f°° f d0id02de3 KiC) tanhC e^^^-R '=°>'h 5 



(27r)3 cosh^ ^ + sinh^ 04 sinh ( 



TT tanh^ fr tanh^ ^^-^ ff — 

2 J--"- 2 -^^ cosh^^ + sinh^6l, 

3 3 



1 



(cosh^ ^ — sinh^ §4 



4 cosh ^ sinh ( 



Y[ (cosh^ C - sinh^ 0,) - 4 cosh^ ^ ^ (cosh^ ^ - sinh^ 0,) Y[ sinh 6lj 

i—l i—1 j^i 

3 3 3 

4 cosh^ ^ J]^ sinh 0i - sinh 6*, J| (cosh^ ^ - sinh^ 0j 



where 0^ and A{uj) are defined by 



04 = Arcoshf ^ - V%osh6l,y A{uj) = ((611,612,613) G K'"^ I V%osh 6*^ < ^ - l|. 



(D3) 



(D4) 



2. Oscillating terms 

-»l.osc , y^l.OSC 



The first correction to ([50)1 is given by €24^^^ + C44°'''^. The usual steps yield for the corresponding spectral function 



= 2anM ^ 



/"^ d6l tanh^ cosh 61 - cosh0 cosh^ ^ - sinh^ 



o2Af_RcoshJ 



27r 7„e, 27r cosh^ ^ + sinh"^ cosh 6* + cosh cosh^ ^ + sinh"^ ^ {uj - M cosh 0^ ~ AP 



Snei^(^ + if) (cosh^ ^ - sinh^ ^) + 2 Jmi4:(^ + i f ) cosh^ sinh^ 



cos(2A/i?sinh( 



(D5) 



ImK + i f ) (cosh^ C - sinh^ 0) - 2mtK{0 + i^) cosh ^ sinh ^ 



sin(2Afi?sinh( 
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where 9 and 9' are defined in ([78]) . We stress that in comparison to (|80|) there appears an extra factor K{£,) e^A-f^^osh^^ 
which shows that (jD5|) is indeed the first correction in K to ([80| . The first correction in K to ([82|l is given by 



2 i°° d( /•"' ^6* i^(Otanh^ g2Affl,cosh? tanh^ ^ 



° io 27r 27r cosh^ ^ + sinh^ 9 cosh^ ^ + sinh^ 9 y/ {uj ~ M cosh 9)^ ~ 
(9\cK{9 + i f ) 5neii:((? + i f ) - 3mK{9 + i f ) 3mii:(^ + i f 

• (cosh^ ^ - sinh^ 6*) (cosh^ ^ - sinh^ ^) - 4 cosh^ ^ sinh 6* sinh 9j 
- 4 fHeK (fi* + i f ) 3mK + i f ) cosh $ (cosh^ ^ - sinh^ 61) sinh 9 + (cosh^ ^ - sinh^ 6*) sinh 

• cos(2Mi?(sinh6' + sinh^)) 

(3mK{9 + i f ) 3mK{9 + i f ) - inei^(6' + i f ) «KeJs:(6' + i f cosh^ 

• ((cosh^ C - sinh^ 9) sinh ^ + (cosh^ ^ - sinh^ 9) sinh 6l) 

- fneif (6* + i f ) 3mK{9 + i f ) (^(cosh^ ^ - sinh^ 9) (cosh^ ^ - sinh^ 9) - 4cosh2 ^ sinh 61 sinh^ 

• sin (2A'fi?(sinh0 + sinh6i)) \. 

(D6) 

The oscillating contribution to local spectral function to second order in K is given by 

STj^iuj, R, r) = S'^rio^, R, r) + Sl^i^. R, r) + R, r) + i?, r). (D7) 

We observe that the integrands are exponentially suppressed for large distances from the boundary. 

3. Contributions from the boundary bound state 

The second-order contribution of the boundary bound state to the spectral function is given by 

S)^^K%oj, R, r) = S',t\^, R, r) + S]:t\^: R, r) + S°t^{uj, R, r) + S°t'{uj, R, r), (D8) 



where 



{Sl^''' + Sli'''') {u, R,r^Q)= Aol cot v tan | e^*^^™^- 



r°° K{£.) tanhCe2*^«™^h« ^^^^^ _ ^^^^ cosh | - sin^; cosh^ ( - sinh^ i (^^) 
Jo 27r ^(cj - i\/sint;)2 - cosh ^ + cos w cosh | + sin u cosh^ ^ + sinh^ ^ 
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and 



2^2 



V 

cot V tan — e 



2MRCOSV 



2 

-R^(C) tanh^e2^^^"°'^'>« cosh^-cosw 



27r ^(w - M sinu)2 - ilP cosh ^ + cos t; cosh^ ^ + sinh^ | 

*Heii:(^" + if ) (cosh 2^'- cos 2w - 2) (cosh^^ - sinh^$")/2 + 4 cosh^ sinh^ cos? 

+ 3mK + i f ) sinh i (cosh 2^ - cos 2v - 2) cosh ^ - 2 cos w (cosh^ i - sinh^ |) 

• cos (2Mi? sinh ^) 

yKzK{i + i f ) sinh I 2 cost; (cosh^ ^ - sinh^ |) - (cosh2| - cos2w - 2) cosh^ 

V 2mK (C + i f ) (cosh 2| - cos 2v - 2) (cosh^ ^ - sinh^ /2 + 4 cosh ^ sinh^ | cos i 

• sin(2Afi?sinhC") I. 



(DIO) 
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